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Abstract. Urban road networks have distinct geometric properties that are partially determined by their 
(quasi-) two-dimensional structure. In this work, we study these properties for 20 of the largest German 
cities. We find that the small-scale geometry of all examined road networks is extremely similar. The object- 
size distributions of road segments and the resulting cellular structures are characterised by heavy tails. As 
a specific feature, a large degree of rectangularity is observed in all networks, with link angle distributions 
approximately described by stretched exponential functions. We present a rigorous statistical analysis of 
the main geometric characteristics and discuss their mutual interrelationships. Our results demonstrate 
the fundamental importance of cost-efficiency constraints for in time evolution of urban road networks. 



PACS. 89.40.Bb Land transportation - 89.75.Fb 
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1 Introduction 

During the last decade, the understanding of the structure 
and dynamics of complex networks and their mutual in- 
terplay has made enormous progress 1 2 3:4 5:6 7 8 19110] , 
A specific type of networks are transportation systems, in 
which the directed flow of individual transportation units 
is often bounded to a spatially discrete substrate consist- 
ing of mutually connected tracks. The merges, diverges, 
and crossings between these tracks can be considered as 
the network nodes, and their physical connections as the 
corresponding links. The resulting graph properties deter- 
mine - together with the link capacities and the actual 
spatio-temporal demand pattern - the efficiency of trans- 
port in the system. 

In modern society, the optimal performance of trans- 
portation systems is of paramount importance for daily 
life. As a consequence, the structure and time evolution of 
such networks often result from careful planning processes. 
Numerous characteristic parameters have been studied for 
rail- and subway [rill2ll3li4ll5li6ll7li8ll9l20l21l22l23l24| . 
airport |25l26l27l28l29l30l3ll32 33 34 35 jg], urban tran- 
sit [37], public transportation |38|39|40|41|42|43|44] , and 
maritime transport networks [45 46 47 on a variety of 
different spatial scales, revealing that many characteris- 
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tics do strongly differ between different systems. Ama- 
ral et al. [25] used this observation for a classification of 
small- world networks (i.e., graphs with a short average 
path length and a high degree of clustering) into scale- 
free, broad-scale and single-scale networks. For example, 
the world airport network has been classified as a single- 
scale network [25 , whereas the node degree distributions 
of the Italian airport network 133] and the Chinese rail- 
way network [3T] are scale-free. Gastner and Newman [?5] 
proposed a connection between the formation of airport 
and highway networks. Although both types of networks 
are intrinsically different, their statistical features can be 
understood as the result of the same optimisation problem 
with one parameter in the goal function varied, leading to 
a variety of network topologies. 

In contrast to many networks in other areas of research 
(e.g., biology, telecommunication, or social sciences), ur- 
ban road systems can be (in good approximation) con- 
sidered as planar networks, i.e., links cannot "cross" each 
other without forming a physical intersection (node) as 
long as there are no tunnels or bridges. Therefore, the 
study of urban road networks allows distinguishing which 
statistical properties are actually specific for planar trans- 
portation systems. Due to reasons of practicality, which 
will be further studied in the course of this paper, and 
given the planarity constraint, it follows that the node 
degrees are practically limited (i.e., the networks cannot 
be scale-free) and the average path length is much larger 
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than for scale-free networks with a comparable number of 
nodes. Both features are therefore distinctively different 
from those of most elementary evolving network models, 
such as (Erdos-Renyi) random graphs |49|50) . (Baraba- 
si- Albert) scale-free [51152153] . or (Watts-Strogatz) small- 
world networks [S] . Even more, the historical development 
of cities can often be understood as the result of cer- 
tain self-organisation processes, i.e., road networks usually 
have not been globally designed following some general 
growth mechanisms of urban areas [54 55 56 57 58 59 60 61], 
which are characterised by scaling laws in various urban 
supply networks and related growth processes [62163)64] . 

While many results obtained so far considered the prop- 
erties of the dual (or connectivity) graphs associated with 
urban road networks (i.e., graphs where roads are repre- 
sented by nodes, which are linked if the physical roads 
do ever mutually intersect |65|66] ). in this contribution, 
we will study and thoroughly compare the characteristic 
geometric properties on the level of the elementary units 
(intersections, road segments, cellular structures) in phys- 
ical space 67J. Unlike in past related studies focussing 
on rather small road networks, the aim of this work is to 
better understand the specific geometric features of the 
complete urban road networks of present-day cities, dis- 
regaring their specific historical development. In order to 
study a homogeneous sample, we consider 20 of the largest 
German cities, which are characterised by rather differ- 
ent traffic conditions and large-scale morphologies (see 
Fig. [I] for some examples). In this paper, we will show, 
however, that the small-scale geometric properties of the 
corresponding road networks show a very large degree of 
similarity suggesting that the common features determin- 
ing the network topology are the planarity and the re- 
quirement of an efficient use of the available space, which 
leads to the preferred emergence of rectangular structures. 

This paper is organised as follows: In Section 2, the 
study area and used data are described. Sections 3 to 
5 summarise our results on the statistical properties of 
nodes, links, and cells as the elementary geometric units 
of a road network. A detailed discussion of our findings 
in relationship to recent results from other studies is pre- 
sented in Section 6. Finally, some outlook on further rel- 
evant directions of research is given. 

2 Study Area 

In this work, we thoroughly extend a recent study |67j on 
the urban road networks of 20 of the largest German cities. 
The large-scale morphology of these networks, which is 
determined by geographical restrictions (hillslope, rivers, 
soil properties, etc.) and the particular historical develop- 
ment (creation of rail- or highways, emergence of certain 
functional requirements, war damages, etc.), differs con- 
siderably between the individual cities (Fig. [I]). 

Our data have been obtained for the year 2005 from the 
geographic data base TeleAtlas MultiNet™ (for detailed 
information, see www.teleatlas.com), which is used in 
numerous applications such as Google™Maps, real-time 
navigation systems (TomTom), and urban planning and 



management. From the 21 largest cities, only Hannover 
has not been taken into account because of a significantly 
lower spatial resolution of the urban area in the initial 
data base, which can be seen from the spatial density of 
intersections (nodes) being by a factor of 5-20 smaller than 
for the other 20 top-ranked cities (see Tab. [I]). Moreover, 
since the initial data base has not been sufficiently homo- 
geneous (e.g., in terms of the representation of separate 
lanes), we have applied a two-step preprocessing: 

In a first stage, multiple and partially directed links 
have been combined by defining that each pair of nodes 
can only be connected by at most one undirected link. Re- 
moving replications and directional information, the num- 
ber of links forming the network (and, hence, the average 
degree of the nodes) is substantially reduced (see Tab. [I]). 
In the following, the resulting undirected network repre- 
sentation will be referred to as the original data. 

In order to better understand some of the geometric 
features of the networks, a second preprocessing step has 
been applied. First, all nodes of degree k n = 2 that repre- 
sent curved or sharply bended roads rather than physical 
intersections of different streets have been removed from 
the data base. This leads to an additional decrease in the 
number of both nodes and links in the networks. Second, 
for convenience, all dead-end roads have been removed 
from the network, since such roads also have only minor 
effects on the traffic conditions in the networks. Formally, 
this has been realised by excluding all nodes of degree 
k n = 1 and their connecting links from the original data. 
Wherever this procedure reduced the degree of other nodes 
to k n = 2, these have also been removed. 

3 Node Properties 

Recent studies on urban road networks have often focussed 
on the statistical properties of the nodes. Apart from their 
spatial density (which is additionally reflected in the link 
and cell properties, which will be discussed in Sections [4] 
and [EJ respectively), nodes are mainly characterised by 
their degree distribution p(k n ), where the degree k n of a 
node in an undirected network is defined as the number of 
road segments (links) that are connected at the respective 
position in geographical space. 

Several authors have already considered the degree 
distribution in different road networks. Buhl et al. |68j 
studied 41 road systems from small- to intermediate-sized 
cities around the world. For the average node degree (k n ), 
they found values between 2.02 and 2.86. The observed 
frequency distributions of the node degrees showed an ex- 
ponential tail for k n > 3 with a variety of possible charac- 
teristic scale parameters, indicating the presence of single- 
scale networks. Scellato et al. [55] compared the proper- 
ties of the road network of the historically grown city of 
Bologna with the planned square- like grid structure of San 
Francisco, yielding average node degrees of (fc„) = 2.71 
and 3.21, respectively. Cardillo et al. [70] argued that this 
observation can be explained by the extraordinarily high 
abundance of "full" intersections (k n — 4) in grid-like 
cities, with p(k n = 4) > p(k n — 3), whereas k n = 3 
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Fig. 1. Large-scale appearance of the road networks of the four largest German cities. The respective scales indicate basic units 
of 1 km length. 



is much more abundant in naturally grown settlements. 
Crucitti and co-workers |71l72l73l74j demonstrated that 
self-organised and planned cities also differ in the proba- 
bility distributions p(h) of the link lengths (unimodal for 
naturally grown, but irregularly shaped for planned cities) 
and the scaling behaviour of different ccntrality measures. 
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Table 2. Absolute and relative frequencies of node degrees k n 
in the urban road networks of German cities, and dependence 
of the mean lengths (li) and associated standard deviations of 
the connected links on the node degrees. 



In the German road networks, we observe that nodes 
with k n — 4 occur much more often than one would expect 



for a single-scale network (i.e., a network with an expo- 
nentially decaying degree distribution, see Tab.^l, which 
is consistent with previous studies. This very high abun- 
dance obviously results from intersections between con- 
tinuing roads. Although the tail of the degree distribution 
resembles an exponential decay, due to the low maximum 
node degree in planar networks, it is hardly possible to 
statistically identify any particular type of scaling |75) . 

We note that the average node degrees (fc„) of the in- 
dividual cities lie within a very narrow range between 3.17 
and 3.31 (see Tab. [3]), although the spatial density of nodes 
differs strongly (Tab. [I]). The average value for all cities, 
(k n ) = 3.25, is remarkably higher than that reported for 
the US interstate highway network ((k n ) = 2.86 [H]), but 
still far below the theoretical limit for planar networks 
((k n ) < 6 [76]). The fact that urban road networks have 
a higher mean node degree than highway networks is pri- 
marily related to the much higher abundance of nodes of 
degree k n = 4. As it will be discussed later, the abun- 
dance of nodes with degrees 3 and 4 is closely related to 
the dominating presence of perpendicular structures in ur- 
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Rank 



2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 



-CitjT 



Population 



Berlin 

Hamburg 

Munich 

Cologne 

Frankfurt 

Dortmund 

Stuttgart 

Essen 

Dusseldorf 

Bremen 

Hannover 

Duisburg 

Leipzig 

Nuremberg 

Dresden 

Bochum 

Wuppertal 

Bielefeld 

Bonn 

Mannheim 

Karlsruhe 



3,392,425 
1,728,806 
1,234,692 
968,639 
643,726 
590,831 
588,477 
585,481 
571,886 
542,987 
517,310 
508,664 
494,795 
493,397 
480,228 
388,869 
363,522 
324,815 
308,921 
308,759 
281,334 



Area 
(km 2 ^ 



753 
311 
405 
249 
281 
208 
210 
218 
318 
204 
233 
293 
187 
328 
146 
168 
259 
141 
145 
173 



37,020 
19,717 
21,393 
14,553 
9,728 
10,326 
10,302 
11,387 
8,237 
10,227 
1,589 
6,300 
9,071 
8,768 
9,643 
6,970 
5,681 
8,259 
6,365 
5,819 
5,505 



n: 



inal 



N, 



87,795 
43,819 
49,521 
29,359 
18,104 
22,579 
21,934 
24,537 
16,773 
21,702 
3,463 
14,333 
21,199 
18,639 
22,307 
15,091 
11,847 
18,280 
13,746 
12,581 
11,866 



50,632 
25,641 
28,579 
18,182 
12,366 
12,512 
12,748 
14,124 
10,377 
12,748 
1,856 
8,141 
11,643 
11,157 
11,977 
8,437 
6,869 
9,880 
7,762 
7,607 
7,004 



13,545 
5,867 
7,074 
3,602 
2,622 
2,176 
2,414 
2,704 
2,114 
2,477 
266 
1,832 
2,561 
2,339 
2,321 
1,455 
1,177 
1,613 
1,381 
1,765 
1,488 



19,931 
9,044 

11,058 
5,395 
3,911 
3,281 
3,612 
4,093 
3,124 
3,827 
411 
2,837 
3,753 
3,543 
3,346 
2,233 
1,750 
2,546 
2,094 
2,674 
2,204 



Processed 



33,000 
14,603 
17,934 
8,724 
6,394 
5,249 
5,848 
6,598 
5,099 
6,106 
640 
4,573 
6,137 
5,762 
5,488 
3,555 
2,803 
4,035 
3,365 
4,364 
3,594 



12,929 
5,458 
6,737 
3,259 
2,445 
1,923 
2,164 
2,444 
1,922 
2,198 
225 
1,710 
2,355 
2,141 
2,109 
1,276 
1,029 
1,461 
1,238 
1,652 
1,366 



pn 

[km~ 21 



"22T4 - 
12.0 
35.6 
13.3 
15.7 
11.7 
17.4 
19.5 
14.3 
12.0 
2.0 
12.2 
12.8 
18.9 
10.2 
15.3 
10.4 
9.8 
14.9 
18.4 
12.7 



Table 1. Characteristics of the urban road networks of the 21 largest German cities in the original and the processed data: 
number of nodes N„, links Ni, and cells N c in the network. The quantity N* corresponds to the content of the original data base 
with partially directed links, whereas Ni reflects the number of links in the reduced undirected network data base. In addition, 
the spatial density p n of intersections in the processed data is given. 
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Wuppertal 
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Table 3. Geometric characteristics of the urban road networks of the 21 largest German cities: node degree k n , link length 
link angle 9i, rectangularity parameters cos 2 29i and p, topological cell degree k c , cell area A c , and form factor (j> c . (•} and a. 
represent mean values and standard deviations taken over all objects in the respective city. Note that the results partially differ 
from those in [67] due to the different preprocessing of the data in our study. 



ban road networks, which is most probably determined by 
both architectural and economic requirements. 

Comparing the results for the German cities with the 
outcomes of other recent studies, it is remarkable that 
both Buhl et al. [68] and Scellato et al. [69] found lower 
values of the mean degree, which can however be easily 
explained by our preprocessing removing all nodes of de- 
gree k n = 1 and 2. If these nodes had been included in 
the data, the mean node degrees would have taken values 
between 2.393 (Bielefeld) and 2.735 (Berlin), which is in 



reasonable agreement with the results of |68l69j . The lat- 
ter values follow directly from Tab.[lj since for undirected 
planar networks, the mean node degree is determined by 

N, 

(U = 2^, (1) 

where Ni and N n are the total numbers of links and nodes, 
respectively. 

The fast decay of the node degree distribution is a 
direct consequence of the planarity of road networks. In 
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terms of practical usage, road intersections having too 
many connections do not provide a feasible solution for 
controlling conflicting traffic flows, fnstead, at such inter- 
sections the total waiting times are significantly enhanced 
and bottlenecks are formed, reducing the overall trans- 
portation efficiency on the network. 

4 Link Properties 

The link-based geometric properties of a road network are 
closely related to the spatial distribution and to the degree 
of the nodes. However, we note that links do not always 
connect neighbouring nodes along their shortest distance, 
which is due to secondary factors like preferred building 
geometry, soil properties, or hillslope. Thus, one has to dis- 
tinguish between the properties of the physical (curved) 
links and those of the Euclidean (shortest-distance) dis- 
tances between their starting and end points. In the lat- 
ter case, all existing links are assumed to be non-curved 
(cf. our preprocessing), resulting in a systematic bias to- 
wards shorter link lengths. In a similar spirit, depending 
on whether the original or processed data are used, the 
term link itself can be identified with either a straight 
road segment without any curves and interior connections 
to other roads (i.e., nodes with k n — 2 exist) or the actual 
physical (possibly curved) connection between two prac- 
tically relevant nodes (k n > 3). In the first case, the links 
are systematically shorter than in the second one. 

Unless it is explicitly stated otherwise, in the following, 
we will consider the properties of the "physical", poten- 
tially curved links. For this purpose, we will define the 
link lengths and directions according to the original data, 
but combining several links into one whenever intermedi- 
ate nodes have been removed in our second preprocessing 
step (i.e., respective node positions and degrees are taken 
from the processed data). 

4.1 Link length distributions 

Travel time is probably the most important indicator for 
the efficiency of a road network. However, this quantity 
strongly depends on secondary parameters such as the 
speed limit and the particular traffic situation (involving 
factors like time of the day, spatial distribution of differ- 
ent functional areas, etc.). Because of this, we consider 
here the travel distance rather than travel time, which is 
measured by the physical length li of each link. 

Let us first discuss the mean link lengths for different 
cities. Table [3] shows that this property varies remark- 
ably (between 128 m (Mannheim) and 199 m (Bielefeld)) 
among the different cities. In particular, there is a ten- 
dency that cities with a small average node degree are 
constructed by somewhat longer road segments. Quanti- 
fying this trend by the linear Pearson correlation coeffi- 
cient r and Spearman's rank-order correlation coefficient 
p as a related distribution-free measure yields values of 
r = —0.55 and p = —0.42, respectively, which supports 
this finding, but underlines that the variations of the mean 




Fig. 2. a) Relative frequency distribution p(h,k n ) = 
p(h\k n )p(kn) of the link lengths of all cities for node degrees of 
k n — 3, 4, 5, and 6 (from top to bottom and black to light grey), 
b) Relative frequency distribution p(h) of the link lengths for 
all 20 cities (black), and upper and lower limits obtained from 
the individual cities (grey lines). 

link lengths are not exclusively determined by the average 
node degrees. 

In order to get a more complete picture, we further- 
more study the full frequency distributions of link lengths 
in the different networks. As a first notable observation re- 
garding the local interrelationships between link and node 
properties, Tab. [2] demonstrates that the link lengths de- 
crease on average with increasing degree of the involved 
nodes. Considering the link length distributions explicitly 
in dependence on k n , Fig. [2] reveals that the qualitative 
behaviour is similar for all node degrees and consists of 
two disjoint regimes: (i) a saturation in the abundance 
of link lengths below 100 m, which seems to result from a 
weakly bimodal shape of the distribution, and (ii) a fat tail 
for longer links. The main differences between the curves 
obtained for different node degrees are a decreasing impor- 
tance of link lengths in the "plateau region" above about 
30 m, and a successive decrease of the maximum value of 
li with increasing k n . Both findings contribute to the ob- 
served decrease of the average link length with increasing 
node degree (Tab. |2j. Our corresponding results confirm 
recent findings of Masucci et al. [77] obtained for the Lon- 
don street network. 

We note that the observed predominance of short links 
coincides with the typical size of building blocks in densely 
urbanised areas (50-100 m). The plateau in the frequency 
distribution is probably related to an inefficiency in having 
very short road sections especially in parts of the city with 
high traffic load, which would result in frequent stops at 
intersections. This hypothesis is supported by the position 
of the (local) maximum of p(h) for k n — 3,4 (possibly 
highlighting the mainly rectangular network structures). 



4.2 Link angle distributions 

In the following, we will discuss the distribution of link 
angles, which encode additional information regarding the 
vectorial properties of road segments in geographical space. 
We note that the meaningful definition of link angles re- 
quires the presence of a planar network, which is assumed 
to be the case in urban road systems. 
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Let £ n (i) denote the set of undirected links that start 
or end at a node i of a planar network. In this case, every 
link j € E n (i) is uniquely defined by three properties: (i) 
the location of node i, (ii) the link length li t j, and (iii) 
the angle ipij with respect to a fixed reference axis (or, 
alternatively, its unit vector). For convenience, this an- 
gle will be counted following the standard mathematical 
convention (i.e., counter-clockwise). For further statistical 
treatment, we will consider the following definition: 

Definition 1: Under the conditions specified above, j r G 
£ n (i) is called right neighbour of j iff 

j r = axgmm{V>i,j' - ipij mod 2w\f € £ n (i)} ■ (2) 

In a similar way, ji € £ n (i) is called left neighbour of j iff 

jl = argmin{^y - ip tjj , mod 2%\j' e £ n (i)} . (3) 

The link angle 6ij is then, for convenience, defined as the 
planar angle between the link j and its right neighbour j r , 
i.e., 

= ^1,3 ~ ^Ijr m ° d 2 7T. ( 4 ) 

The link angle directly relates the spatial distribution 
of the nodes (in particular, their degrees) with the sta- 
tistical features of links, which is a specific characteristic 
of planar networks. Unlike the link lengths, link angles 
exclusively refer to the respective node since 



E < 

J€£„(i) 



•id 



2tt. 



(5) 



As for the link lengths, all 20 road networks consid- 
ered in this study show again approximately the same fre- 
quency distribution of the link angles as displayed in Fig. 
[3^i. This universal behaviour supports the hypothesis of a 
common construction principle of all studied urban road 
networks. Indeed, the mean link angles for the 20 indi- 
vidual cities scatter by far less than the link lengths (i.e., 
between 108.6° (Berlin) and 113.1° (Bremen), see Tab. [3}. 
Note that there is an almost perfect (negative) statistical 
dependence between the mean link angle and the mean 



node degree (fe n ), since 



360°/fc„ 



In general, the link angle distributions p(9i) have two 
pronounced peaks at 90° and 180°, where link angles of 
about 90° are more abundant and show a larger disper- 
sion than those of about 180°. This observation is directly 
related to the node degree distribution discussed above: 
When considering only nodes of degree k n — 4, there is 
only one sharp peak at 90° (see Fig. [3(3), which corre- 
sponds to almost perpendicular intersections of pairs of 
roads as an extraordinarily probable way of forming nodes. 
For fc„ = 3, the link angles around 90° occur twice as of- 
ten as those around 180° (note that the ratio between 
the corresponding peak frequencies in Fig. [3]d is clearly 
smaller than 2, since the dispersion around the 90° peak 
is larger), which implies that the corresponding nodes are 
mainly formed by a straight primary road from which a 
secondary one splits perpendicularly. 




90 180 270 360 

e,[°] 



90 180 270 360 
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Fig. 3. a) Average frequency distribution p(0i) of the link 
angles from the processed data taken over all cities (black), 
and upper and lower limits obtained from the individual cities 
(grey), b) Conditional distributions p(9i\k n ) for node degrees 
k n — 3, 4 and 5 (from black to light grey). c,d) Same as a,b) for 
the double angle distribution p(Si) (only k n = 3,4 are shown 
in d). e,f) Distributions of the parameters cos 2 28i (e) and fi 
(f ) measuring the local deviation from a perfectly rectangular 
profile. The black line again represents the mean taken over all 
cities, grey lines indicate the minimum and maximum values 
for the individual road networks. 



In order to better understand the observed behaviour, 
remember that in the construction of urban road networks, 
minimisation of construction and maintenance costs as 
well as maximum usage of the road have to be taken into 
account [48] . In terms of geometry, the shortest Euclidean 
connection is the best solution for constructing the new 
road at minimum costs. This shortest line drawn from a 
node to an existing link naturally forms a perpendicular 
intersection. If this new road is further continued, four 
right angles are formed. If it touches only, two right an- 
gles and one 180° angle emerge. Moreover, in terms of 
practical usage of the available space (in particular, re- 
garding the typical building geometry) , rectangular shapes 
are preferable. Finally, also in terms of turning vehicular 
traffic, symmetric 90° angles are a particularly reasonable 
choice, since for 9i 56 90°, turning processes become much 
more complicated in terms of curve radius and speed. The 
link angle distribution strongly support this explanation, 
which identifies a main driving mechanism for growing ur- 
ban road networks. 

For a quantitative evaluation of the rectangularity of a 
given road network, we consider some statistics that mea- 
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sure the local deviation from a perfectly rectangular grid. 
As two properly normalised parameters, the squared co- 
sine of the doubled link angles cos 2 29i and the normalised 
minimum absolute deviation from right angles 



a = — mm 
7T fcez 



k 

— 7T 

2 



(6) 



have been computed for all link angles. The mean values 
and associated standard deviations for all cities are listed 
in Tab. [3] According to these values, Munich, Mannheim 
and Berlin have the "most rectangular" road networks. In 
turn, the road networks of Wuppertal and Bochum show 
the least pronounced perpendicular splitting of roads. How- 
ever, all values still lie within a relatively small range. As 
expected, both measures strongly correlate (r = —0.981, 
p = —0.985 for the respective mean values of the indi- 
vidual cities (in a perfectly rectangular network, we have 
(cos 2 26[) = 1 and fi = by definition). However, if both 
characteristics are computed for all nodes in the studied 
networks, we observe a significant deviation from the ideal 
case, which results in a broadening of the corresponding 
probability distribution functions (Fig. (3fe and f). 



4.3 Double-angle distributions 

In the link angle distributions, we have still observed de- 
viations from the perfectly rectangular alignment of inter- 
secting roads, which confirms daily life experience. Hence, 
our results obtained so far do not necessarily imply that 
intersections are also formed by straight road segments. 
In order to further study the straightness feature explic- 
itly, the joint distributions of neighbouring link angles are 
studied next. This analysis can be simplified by introduc- 
ing the following auxiliary variable: 

Definition 2: Under the conditions of Definition 1, the 
double angle Sij associated with a link j that enters a 
node of degree k n > 2 is defined as the angle between its 
left and right neighbours: 



l,3r 



mod 2tt. 



(7) 



While the distribution of the individual link angles 0i 
quantifies the rectangularity of intersections, at least for 
nodes with degree k n = 3 and 4 (i.e., about 99.7% of all 
intersections, see Tab. [2]), the associated double angles Si 
measure the straightness of the crossing roads. It is trivial 
to demonstrate that 



(8) 



The distributions of double- angles (Fig.[3j: and d) show 
the behaviour to be expected from the link angles. In par- 
ticular, for nodes with degree k n — 4, there is only one 
sharp peak around 180°, while for nodes with k n = 3, 
there are two peaks at 180° and 270°. Note, however, that 
although one would expect double angles of about 270° to 
occur twice as often as those of about 180° for k n = 3, we 
observe that the corresponding peak frequencies differ less 



from each other. This again implies that the dispersion is 
larger around the 270° maximum (see Fig. [3]i). Note that 
contributions due to nodes with degrees k n > 5 do not 
have any remarkable impact on the full double-angle dis- 
tribution due to their very low abundance. 



4.4 Statistical model for link angle and double-angle 
distributions 

In the following, our above considerations will be supple- 
mented by a simple statistical model for the observed link 
angle and double-angle distributions. In particular, for 
k n = 3 and 4, the sharp and partially symmetric peaked 
structure of the observed frequency distributions suggests 
that a stretched exponential (or generalised Laplacian) 



P(0i\a,r) 



1 

N 



exp 



(9) 



(or a similar expression for Si) could provide a reasonable 
statistical model. In the latter expression, 

2r f°° 
N=—r{\/a) with r(t)= / e~V -1 tfu (10) 
a Jo 

is a normalisation constant, 9i$ the peak location, r mea- 
sures the dispersion of the distribution around the peak, 
and the exponent a mainly determines the scale of the 
probability decay as 6i deviates from the peak position. 

Stretched exponential distributions have been recently 
used for describing the cumulative frequency distributions 
of a variety of natural as well as economic time series 
[78] (e.g., peak-trough patterns of commodity prices [79]). 
There are only few examples for applications as models for 
probability densities, e.g., of relative price increments (as- 
set returns) in financial time series [80] . In a wider sense, 
generalised Laplacians have been used for approximating 
the distributions of wavelet coefficients in image process- 
ing applications [81] , 

The full distributions of both link and double angles 
show minor asymmetries around their respective peaks, 
which is related to the presence of nodes of degree k n > 4. 
However, for k n = 3, 4, one could expect a symmetric 
behaviour. Note that the full distributions p(0i) and p(5i) 
can be expressed as finite mixture models [82 with compo- 
nents corresponding to different fc„, the statistical weights 
7r^ fen -' [k n =3,4,...) of which are determined by the rel- 
ative frequencies of the different node degrees p(k n ), i.e., 
7r(fc») =p(k n )- k n /(k n ). 

For k n = 4, the double angle distributions for S t < 180° 
and Si > 180° pass different statistical tests for homogene- 
ity (Kolmogorov-Smirnov, Kuiper, % 2 ) on a 99% confi- 
dence level. For the corresponding link angle distributions, 
tests on the explicit data (Kolmogorov-Smirnov, Kuiper), 
however, reject the null hypothesis of symmetric behaviour 
around 9i = 90°. In contrast, since rounding errors in the 
data base could matter here, binning of the corresponding 
data seems reasonable, which leads (for a binwidth of 0.5°) 
to x 2 values that do not allow rejecting the null hypothe- 
sis of symmetric tails anymore. Figure [^demonstrates this 
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Fig. 4. a) Empirical frequency distributions p(\0i — 0i,o\) 
(solid lines) for k„ = 4 (6>; i0 = 90°) at 6 l < 6 lfi (black) 
and 6i > #;,o (grey), and corresponding stretched exponen- 
tial fit (dashed lines), b) Same as a) for the double-angles Si 
{8i,o = 180°). c,d) Same as a,b) for k n = 3. e,f) Same as c,d) for 
9i,o = 180° and J;,o = 270°, respectively. Note that some of the 
lines are hardly visible, since they coincide almost completely 
with others. 



symmetry and underlines that at least for \9i — 90° | > 2° 
and 1 5i — 180° | > 2°, stretched exponentials fit the be- 
haviour of the empirical distribution functions reasonably 
well. Note that the explicit estimation of the correspond- 
ing model parameters is a challenging task, which is be- 
yond the scope of the presented study. However, as a gen- 
eral result, we summarise that both the associated scaling 
exponents a and the dispersion parameters r are clearly 
larger for the link angles than for the double angles. This 
implies a broader distribution near the peak, but a faster 
decay in the tails of p(6i) (see Fig. [31. 

In contrast to the previously considered case, for k n — 
3 the link and double-angle distributions have two peaks 
with different properties. Although the distributions of the 
tails associated with these peaks do not pass the homo- 
geneity tests mentioned above, it can be seen from Fig. [4] 
that there is a considerable degree of symmetry around the 
90° peak oip{9{) and the 270° peak of p(Si), respectively, 
whereas the tails associated with the 180° peaks show 
strong asymmetries. Similar statements hold for the possi- 
bility of describing the empirical distributions by stretched 
exponential functions. The observed asymmetry around 
the 180° peak is possibly related to the mutual overlap 
of the tails associated with different peaks, which cannot 
be sufficiently resolved by the heuristic cutoffs considered 



here in the middle between two peaks. A more sophisti- 
cated statistical treatment might help improving the cor- 
responding results in future studies. 



4.5 Radial distributions of link properties 

The decrease of the typical link lengths with increasing de- 
gree of the involved nodes suggests a possible relationship 
with the spatial distribution of nodes and links within the 
studied cities. Even though a detailed investigation of the 
historical development of the 20 networks is beyond the 
scope of this study, in the following, we will approximate 
the temporal evolution by the distance of an intersection 
from the geometric city centre x. Although this approxi- 
mation is only reasonable for symmetrically grown cities, 
we argue that historical quarters are typically located close 
to the geometric city centres, whereas modern parts can 
more often be found in the outer regions. 

We speculate that nodes with high degrees are usually 
located in densely urbanised areas, in particular, more of- 
ten relatively close to the city centres. In a similar spirit, 
we also expect shorter streets being typical for historical 
city centres. In order to examine this hypothesis, for every 
city, the distances of all nodes from the respective geomet- 
ric centre have been computed. Fig. [5^, demonstrates that 
there is indeed a tendency that high-degree nodes can be 
found more often (in terms of relative frequencies) close to 
the city centres than such with lower degrees. This result 
is again in good agreement with recent findings by Ma- 
succi et al. for the London street network, who observed 
an approximately linear decrease of the mean node degree 
with an increasing distance from the city centre [77]. In a 
similar way, the quantiles of the probability distributions 
of the link lengths show a clear tendency towards shorter 
links near the city centres (Fig. ^jp) , which also supports 
the results for the London street network [77]. However, 
while longer road segments occur especially often at inter- 
mediate distances (about 5-10 km) from the city centres, 
in the outer parts of the examined cities we again find 
a tendency towards a lower frequency of very long links, 
which could be related to a predominance of residential 
areas in these regions, which are typically characterised 
by short roads. 

In addition, we evaluate a possible dependence of the 
road geometry on the distance from the city centres. As 
a parameter describing the curvature of individual road 
segments, we consider the ratio between the corresponding 
link length li and the Euclidean distance between the two 
connected nodes, 



9i 



X; i - X;.2 



> 1. 



(11) 



Figure [HJ: shows that there is actually a tendency towards 
only few links with significant curvature (g\ > 1) close to 
the city centres, which implies that in the corresponding 
parts of the cities, both short and straight road segments 
are particularly common. In contrast, there seems to be no 
clear trend towards more rectangular structures in these 
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Fig. 5. a) Conditional cumulative frequency distributions 
P{r n \k n ) of node distances from the estimated city centres for 
k n = 3, 4, 5, and 6 (from black to light grey). b)-d) 95% to 5% 
(in steps of 5%, from top to bottom) probability levels of the 
empirical conditional distribution functions of the link lengths 
(b), the curvature parameter gi (c), and the rect angularity pa- 
rameter (d) in dependence on the distance r„ from the city 
centres. 

areas as shown by the corresponding distribution of the 
parameter \x in dependence on r n (Fig.[5jl). 



5 Cell Properties 

Roads divide urban areas into cellular substructures. The 
two-dimensional objects, which are formed by closed loops 
consisting of different links, are thus called the cells of the 
road network [57] . Formally, let V c (i) be the set of nodes 
that belong to a given cell i. The set £ c (i) of links forming 
the boundary of this cell is then given as 

Sc(%)= IJ ( U n£„(ro')) J • 

m£V c (i) \m' £V c (i) ,m'^m / 

There is a variety of characteristic geometric properties 
describing a cell, in particular 

— the cell area A c , 

— the topological cell degree k c (i.e., the number of neigh- 
bouring cells [83] , where neighbours are defined accord- 
ing to the presence of a common edge) , 

— the geometric cell degree n c (i.e., the number of straight 
road segments forming the cell, note that Masucci et al. 
[77] called this quantity the cycle length CI), 

— the (maximum) cell diameter d c , and 

— the cell perimeter |84l85j 




Fig. 6. Relative frequency distributions of the topological cell 
degrees k c (a), geometric cell degrees k c (b), cell areas A c (c), 
and form factors <p c (d). The black line represents the mean 
taken over all cities, and grey lines indicate the minimum and 
maximum frequencies obtained in the individual networks for 
given k c , k c , A c and tj> c , respectively. 

It is obvious that the different measures are not mutu- 
ally independent. For example, in order to characterise 
the shape of a cell, it is therefore convenient to combine 
cell area and diameter into one normalised parameter, the 
form factor (or structure factor) <fi c [67 84 85 , which is de- 
fined as the ratio between the actual cell area A c and the 
area of the smallest possible circumscribed circle. Since 
the latter one is determined by the maximum cell diame- 
ter as 

nd 2 

A mm (dc) = (14) 
the defining equation of the form factor reads 

A c AA C 

A m in{d c ) TTC^c 

In the following, the distributions of cell degrees, cell 
areas, and form factors will be studied as examples of the 
typical properties of urban road networks. Since we are 
interested in the physical shapes of the cells, the original 
data will be used here again instead of the processed ones. 

5.1 Cell degree distributions 

In contrast to the nodes, the possible degrees of the cells 
(in both topological and geometric sense) are not restricted 
to very low numbers. Figure^ shows that among the 20 
German cities, there is again a common behaviour of the 
topological cell degree distribution p(k c ) with a fat tail for 
k c > 4. Similar observations can be made for the associ- 
ated geometric cell degree n c (n c > k c , see Fig. [sJd) . The 
average topological cell degrees are of the order of 5, with 
only rather weak variations between values of 4.95 (Karl- 
sruhe) and 5.34 (Bielefeld) as shown in Tab. [3] The mean- 
ing of the topological cell degree can also be illustrated 
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by considering an adjoint graph representation based on 
these two-dimensional objects, i.e., a network where the 
nodes represent the cells of the physical road networks and 
the links stand for common boundaries between the cells. 
It follows from the above results that such adjoint graphs 
form approximately scale-free networks, since the topolog- 
ical cell degree k c in the physical space corresponds to the 
node degree £;* in the associated adjoint representation of 
the cells. 



that among the cells with large degree, there is a signif- 
icant subset characterised by strongly elongated shapes. 
Such cells may be especially observed along geographical 
boundaries (e.g., rivers, highways, or railways) or in spe- 
cific functional areas (e.g., large parks or industrial areas). 
Our aforementioned findings are supplemented by rather 
weak negative correlations between the average topologi- 
cal cell degrees and form factors of the 20 considered road 
networks (r = —0.31, p = —0.31). 



5.2 Form factor distributions 



5.3 Cell area distributions 



Concerning the form factors, we again find very similar 
distributions for all cities (Fig. [6jl) , with preferred values 
between about 0.3 and 0.6. These numbers are in good 
agreement with the maximum possible values for cells with 
degree k c = 3 and 4, 



'(Jfc c = 3) = ^U IUI. 



-in 



0.64, 



(16) 
(17) 



where the latter one is the most abundant type (Fig. [6k) . 
Among the different cities, the average form factor of all 
cells varies (see Tab.[3| in a range between 0.352 (Bremen) 
and 0.411 (Dresden and Mannheim), and the standard 
deviation of <p c for the individual networks takes values 
between 0.159 and 0.175. 

Considering the interrelationships between the differ- 
ent cell characteristics, we expect a dependence of the form 
factor <f> c on the cell degree k c . In particular, </> c is bounded 
from above by the maximum area of a cell formed by a 
given number of nodes. For a prescribed perimeter p c , this 
maximum is taken in the case of a regular fc c -polygon. 
Thus, we find 



kc . ( k c 
— sin — — — 

TT \ 2k r 



1 



(18) 



with <fi" lax — > 1 for large k c . For the German cities, there 
is, however, hardly any cell with <fi c > 0.64 (i.e., the value 
of <fi™ ax for k c = 4). This finding underlines that even cells 
with a high degree mainly have rectangular shape, which 
is again a consequence of the preferred form of buildings 
and specific functional areas within urban settlements also 
reflected in the distributions of link and double-angles. 
Even more, Fig. shows that the most frequent form 
factors are even smaller than 0.64, which suggests that 
many of these rectangular cells of high degree have elon- 
gated rather than squarish shapes. 

In order to further support the latter interpretation, 
the dependence of the form factor distribution p((p c ) on the 
topological cell degrees has been explicitly studied. From 
the results shown in Fig. [7^,, we identify an increasing 
probability of very low form factors for high cell degrees, 
i.e., cells with a possibly more complex (less rectangular) 
shape. Moreover, there is hardly any cell with a form fac- 
tor above 0.64 even for k c 3> 4. Both observations suggest 



Like for all other quantities studied so far, the distribu- 
tions of cell areas also show a common behaviour for all 
considered German cities (see Fig. (6k). In particular, the 
qualitative characteristics are similar to those of the link 
lengths, with almost constant cell size probabilities below 
about 10 4 m 2 and a fat tail for larger cell areas. However, 
the behaviour at small sizes is distinctively different, i.e., 
small cells are observed more frequently than those of in- 
termediate size, whereas in the link length distributions, 
a plateau-like saturation behaviour has been found (see 
Sec. One reason for this is the influence of the cell 
shape: there are many cells that are exclusively formed by 
short road segments, some arising from a combination of 
short and long links (i.e., with elongated structures), and 
only a few cells exclusively resulting from very long links. 
In the following, the corresponding interrelationships will 
be discussed in some detail. 

First, we note that for a known link length distribu- 
tion, the cell areas are closely related to the cell degrees, 
the corresponding relationship being mediated through 
the associated form factors. For the respective mean val- 
ues computed for the 20 individual cities, we actually find 
that the relationship between cell areas and cell degrees is 
characterised by a linear correlation coefficient of r = 0.76 
(p = 0.67), which is remarkably high. 

In order to better understand the dependence of the 
cell area on the cell degree, let us assume that all links 
forming the network would have the same length L (of 
course, this is only possible for rather generic network 
structures such as triangular, square, or hexagonal grids). 
In such a case, the maximum possible cell area (which is 
occupied by a regular A: c -polygon) is a function of the cell 
degree k c as 



A, 



L 2 

;{kc) = ~^~^ c 



1 



tan 



(*-) 



(19) 



and approaches asymptotically the corresponding expres- 
sion for the area of the circumscribed circle, 



4 • 2 

* sin 



4 TT ' 



(20) 



From the latter expression, it follows that the rela- 
tionship between cell areas and link lengths is at most 
quadratic. However, in real-world urban road networks, 
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Fig. 7. a,b) Frequency distributions (colour-coded) of form factors <j> c (a) and cell areas A c (b, in logarithmic scale) in dependence 
on the topological cell degrees k c . The white solid line indicates the theoretical maximum value according to Eq. ( 18 1, while the 



black lines correspond to the mean values as a function of the cell degree. c,d) Frequency distributions (colour-coded) of the 
form factors in dependence on the cell area (in logarithmic scale) for all cells (c) and for cells with fc c = 4 only (d). Black lines 
again indicate mean form factors in dependence on the cell area. 



such a dependence cannot be found. On the one hand, 
even simple rectangular cells of degree k c — 4 show a non- 
trivial distribution of the aspect ratio A c (which is related 
to the form factor as cj) c = 47r~ 1 A c j '(Z\ 2 + 1)). In partic- 
ular, because the assumption of constant link lengths is 
violated, systematic deviations from the quadratic depen- 
dence arise. On the other hand, since even cells of higher 
degree typically have rectangular rather than higher-order 
polygonial shape, the typical cell areas scale clearly sub- 
quadratically with k c (Fig. \7p). An alternative explana- 
tion for this observation could be that cells with a large 
degree have a tendency towards being formed by links 
with shorter lengths. However, the average values for the 
20 German cities indicate that there is in contrast a sig- 
nificant positive correlation between link lengths and cell 
degrees (r = 0.70, p — 0.52) as well as cell areas (r = 0.93, 
,9 = 0.83). 

Finally, we observe a non-trivial dependence of the 
form factors on the corresponding cell areas, which is shown 



in Fig. [7j:. In particular, it turns out that cells cover- 
ing a large area have a tendency towards more squarish 
shapes (i.e., larger form factors), while for smaller cells 
(A c < 10 4 m 2 ), the average form factors decreases signifi- 
cantly, which implies more elongated structures. This in- 
terpretation is supported when considering only the subset 
of cells with a topological cell degree of k c = 4 (Fig. [7]i). 
However, from the latter subset, we find that very small 
cells (A c < 10 3 m 2 ) again show a clear tendency towards 
squarish shapes (i.e., higher <p c ). These findings show that 
the typical aspect ratio of cells in urban road networks 
is related with the cell size in a complicated way. Follow- 
ing our previous considerations, we note that the empiri- 
cally observed dependence is most probably determined by 
both building sizes and transport efficiency. In particular, 
we argue that the transport efficiency decreases for small 
cells (due to frequent stops of traffic) as well as for large 
cells (increase of travel distance due to the predominant 
rectangular structures). In this respect, the emergence of 
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r c [km] r c [km] 



Fig. 8. 95% to 5% (in steps of 5% from top to bottom) prob- 
ability quantiles of topological cell degree k c (a), geometric 
cell degree k c (b), cell area A c (c), and cell perimeter p c (d) 
obtained from the 20 considered cities in dependence of the 
distance of the centroid midpoints r c from the city centres. 

more elongated shapes for cells of intermediate size could 
be regarded as a certain trade-off, including parallel roads 
with a small spacing and perpendicular roads with larger 
mutual distances as a typical spatial pattern. The fact 
that small cells approach again less elongated shapes is 
probably related to the fact that these are often located 
in certain functional areas (e.g., residential areas) with 
a rather low traffic volume, which would imply that the 
dependence of cell shapes on cell sizes is at least partly 
influenced by the expected traffic conditions. 



5.4 Radial distributions of cell properties 

In addition to the aforementioned relationships between 
cell and link properties on the one hand, and link and 
node properties (cf. Sec. [4]) on the other hand, we note 
that it is possible to also evaluate correlations between the 
node degrees and the areas and degrees of the surrounding 
cells. From the average values of the 20 German cities, 
we find r = -0.53 (p = -0.44) and r = -0.83 (p = 
—0.86), respectively. The last result is a consequence of the 
fact that cells of large degree have most often rectangular 
shapes and are formed by sequences of links that connect 
nodes with degrees of 3 and 4. 

Since both node and link properties depend on the 
spatial distance from the city centres (cf. Sec. H), we ex- 
pect a corresponding dependence for the cell properties 
as well. Masucci et al. [77] recently found that for the 
London street network, both average values and standard 
deviations of cell degree and cell area increase towards the 
outer parts of the city. Qualitatively similar findings have 
been reported for the distributions of areas and perimeters 
of individual buildings [55] (which are closely related to 
those of the cells of the associated urban road networks) . 
The aforementioned results coincide with our observations 



made for the 20 German cities, which can be seen in Fig.[H] 
However, in contrast to the London case, we observe (as 
for the link lengths) a saturation of cell degrees, areas, and 
perimeters at intermediate distances from the city centres 
of about 5-10 km. This difference could be related to a 
possibly different structure of the analysed data as well as 
actual differences in the urban patterns. 

As a final remark, we emphasise that the spatial het- 
erogeneity of cell properties is related to the non-trivial 
large-scale morphology of the networks (Fig. [I]), which is 
reflected in the differences of global characteristics like 
the population density or the spatial density of intersec- 
tions p n (see Tab.[T]), the latter of which is directly linked 
with network properties like mean link length or mean 
cell area. Among the studied networks, we particularly 
find as the two extreme cases that (i) Bielefeld has many 
large cells of high degrees, but a low average node degree 
(i.e., many k n — 3 nodes), whereas (ii) Mannheim has 
much fewer large cells and, thus, a lower mean cell de- 
gree. The different relative frequencies of large cells most 
probably indicate the presence of distinct functional areas 
(industry, business, parks, etc.) or geographic constraints 
(mountains, lakes, or rivers, etc.). 

6 Discussion 

Recent investigations on the properties of road networks 
in physical space have often been restricted to node prop- 
erties (i.e., distributions of degrees and other centrality 
measures) discussed for small parts of urban areas only 
|68I69I70I71I72I73I74| . These studies revealed distinct dif- 
ferences between the network properties of historically 
grown (self-organised) and planned cities, as discussed in 
Sec. 3. In contrast, we have studied various geometric 
properties of the complete urban road networks of 20 Ger- 
man cities, which strongly differ among each other in terms 
of their historical development and large-scale morphol- 
ogy. Nevertheless, a large degree of universality is found 
in all considered properties, which reflects common de- 
sign principles in both historic and modern cities such as 
maximum area usage, optimisation of traffic conditions 
(with respect to necessary stops and turning processes), 
and minimum construction and maintenance costs. The 
combination of the aforementioned principles implies the 
presence of perpendicular structures [48] as supported by 
our empirical findings. 

In the last years, studies on the network properties 
of urban road systems in physical space have been supple- 
mented by a number of investigations addressing the topo- 
logical features of the corresponding connectivity graphs. 
For example, Jiang and Calamund }65j found that urban 
street networks can be classified as small-world networks 
with small average path length and high degree of cluster- 
ing. Kalapala et al. |75j reported a universal scaling be- 
haviour p(k n ) ~ fc~ 7 in the degree distribution p(k n ) ob- 
tained from the dual representations of the national road 
networks of the United States, England, and Denmark, 
with exponents 7 of 2.2 to 2.4. This result is in certain 
contradiction to findings of Jiang and Calamund [65] and 
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Rosvall et al. [87], who did not observe power-laws in the 
connectivity distributions obtained from small-scale ur- 
ban road networks. Wagner [5S] found that the distance 
distribution in the dual graphs is described by a rescaled 
Poissonian, the parameters of which are closely correlated 
with the scaling of the power-law distribution of the rank- 
ordered street lengths. 

The contradicting results on the emergence of scal- 
ing in the connectivity distributions of the dual graphs 
can again (at least in some parts) be understood by a 
distinction between self-organised and planned cities |70j . 
Porta et al. [71] showed that a power-law with an exponent 
of 7 = 2.5 ±0.1 describes the degree distribution of the 
road network of the naturally grown city of Ahmedabad, 
whereas other networks studied in their work showed ei- 
ther no scaling at all or included too few intersections to 
conclude with a sufficient significance about the presence 
of a power law. Jiang [89 found a distinct scaling in a set 
of 46 urban road networks with scaling exponents 7 be- 
tween 1.7 and 2.6. In addition, he reported that in planned 
cities, very long streets (i.e., nodes of the dual graph with 
a high connectivity) result in a second scaling region with 
additional exponents of about 3.4 to 4.2. 

From the aforementioned results, we conclude that the 
fact that historical urban structures have grown over much 
longer time intervals makes their statistical features being 
rather specific. The set of German cities studied in this 
contribution can still be considered as a relatively homo- 
geneous sample in the sense that the geographical condi- 
tions of the study area, the total age of the cities, and the 
economic and demographic history determining the evo- 
lution of the urban areas are sufficiently similar. This is 
also reflected by the universality of the geometric proper- 
ties reported in this paper. In turn, this argument implies 
that the presented results do not provide enough empir- 
ical arguments for an emergence of comparable statisti- 
cal properties in cities with completely different historical 
or cultural background. It will therefore be an important 
question for future research to which extent the presented 
results can be generalised independent of the region or his- 
tory of a city. Moreover, for the German cities studied in 
this work, we can hardly conclude whether or not histori- 
cal and modern parts have comparable geometric proper- 
ties, since there is no clear distinction between both types 
of areas due to multiple small historical centres (e.g., in 
the case of Berlin) , past damages of historical areas due to 
wars and governmental policies, and other related effects. 

From a conceptual point of view, the structure of nat- 
urally grown urban systems has been proposed being the 
result of some generic self-organisation processes similar to 
those observed in biological systems. For example, Schaur 
|54j presented a comparison of some simple network prop- 
erties of human settlements with those of biological sys- 
tems such as dragonfly wings or maple leaves. Buhl et al. 
|90j performed extensive experimental studies on the de- 
velopment of ant galleries in restricted planar areas. They 
found a variety of different network patterns with average 
node degrees (k n ) between 2.0 and 3.06 and single-scale 
exponential degree distributions. Similar results have re- 



cently been reported by Perna et al. |91j for three-dimen- 
sional gallery networks of termite nests (with (k n ) between 
2.06 and 3.27). Due to our preprocessing for reducing the 
road networks to their geographically relevant structures, 
we are not able to systematically compare the results for 
the German cities with those obtained for other (biolog- 
ical as well as technological) transportation networks. A 
corresponding analysis will be subject of future studies. 

We emphasise that in the framework of the presented 
purely geometric analysis, we are not able to definitely 
distinguish between small-scale structures resulting from 
either historical (self-organised) growth or careful urban 
planning. Since both factors can be expected to apply 
in our case study, our findings suggest that some funda- 
mental local geometric properties of urban road networks 
are indeed almost independent of the particular history 
of a city. While the relevant spatial dimensions (i.e., link 
lengths and cell areas) of the basic units show a substan- 
tial scatter between the studied road networks, the corre- 
sponding distribution functions are - qualitatively - ex- 
tremely similar. Therefore, we speculate that the phys- 
ical sizes of basic units can show some dependence on 
the historical evolution of the network (e.g., smaller cells 
and shorter links in historical city centres as revealed by 
our analysis). In a similar way, specific geographic con- 
straints and other secondary factors matter here, whereas 
fundamental principles related with cost and transporta- 
tion efficiency apply to both historically grown as well 
as strategically planned cities. Accordingly, the identified 
"universal" geometric properties do not vary significantly 
among the road networks under study, although the par- 
ticular history of the corresponding urban structures may 
be quite different. 

Despite the limitations discussed above, the large de- 
gree of similarity in the geometric properties of urban road 
networks found in this study clearly supports the assump- 
tion of some universal design principles beyond the evo- 
lution of urban structures, which seem to apply to both 
historically grown as well as planned road networks. Since 
road capacities have not been considered in our investiga- 
tions, one could assume a certain independence of struc- 
tural properties from the traffic volume based on these 
first principles. However, we have observed certain non- 
trivial interrelationships between different geometric prop- 
erties (reflected by the complex individual large-scale mor- 
phologies of different cities), which suggest that traffic con- 
ditions actually matter in the growth of urban structures. 
Hence, in the spirit of recent results regarding the self- 
organisation of urban settlements, we argue that there is 
a co-evolution of structure and dynamics (even beyond the 
transportation capacity of roads) already on the geomet- 
ric level, which is supported by daily life experience. This 
implies that urban road networks have to be considered as 
adaptive co-evolutionary networks (92193j . the geometric 
properties of which arise from both self-organisation prin- 
ciples and planning strategies. We conclude that there is 
a considerable importance of first principles such as min- 
imum travel times and maximum area usage in the evo- 
lution of urban structures, which lead to some universal 
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structural properties of the resulting networks mainly de- 
termined by the planarity constraint [77] . These general 
results, however, need to be further validated by simula- 
tion studies of corresponding sophisticated mathematical 
models of growing networks. 

One particularly interesting recent result on possible 
relationships between road network geometry and traffic 
dynamics has been given by Lammer et al. |67j , where the 
authors reported estimates of effective travel time-based 
dimensions of d > 2 for Germany's largest cities. We note 
that this finding cannot be directly related with the out- 
come of our geometric analysis, since the definition of the 
dimension considered in [57] has been based on traffic dy- 
namics instead of geometric network structures. Concep- 
tually, we suggest that higher "dynamic dimensions" could 
be related to a strong heterogeneity in the distribution of 
travel times, which can be triggered by bottlenecks such 
as few bridges connecting different parts of a city (that 
have high internal connectivity), i.e., the presence of well- 
expressed modular structures in the network. In this spirit, 
the dynamic dimensions of |67j are related to meso-scale 
network properties rather than local geometric features of 
urban road systems. 

7 Outlook 

We have demonstrated that urban road networks are a 
class of transportation systems with distinct geometric 
properties (in geographical space), which are controlled 
by planarity as well as cost and transportation efficiency. 
The observed similarity of the distributions of various ge- 
ometric quantities among a set of cities with heteroge- 
neous large-scale morphology suggests that some of these 
features, in particular, the predominance of rectangular 
structures and the emergence of fat tails in object-size dis- 
tributions, result from universal construction principles, 
which calls for additional modelling studies on this phe- 
nomenon. In turn, recent results indicated that at least 
the presence of fat tails in the cell area distributions could 
arise as a natural consequence of the planarity of urban 
road networks [77] , since similar properties can be recov- 
ered for rather general planar network models. Hence, 
it will be an important question for further (both an- 
alytical as well as modelling) studies to distinguish be- 
tween properties directly resulting from planarity con- 
straints and those that are originated in dynamically rel- 
evant self-organisation mechanisms beyond the network 
growth. In this respect, there is also a need for deeper 
studies on potential relationships between the emergence 
of allometric scaling, which is related to power-laws in the 
rank-order distributions of geometric objects [86], and the 
fat tails of physical object-size distributions observed in 
both empirical data and simple models of urban growth 
|67|77|84(85] . In particular, it has to be studied in some 
detail whether or not the corresponding tails actually cor- 
respond to power-laws with a universal scaling exponent 
of 2 for (at least a certain class of) planar networks, as it 
has been recently suggested [77] , and to which extent this 
is reflected by the corresponding link length distributions. 



As a second main line of future research, we empha- 
sise the need of a spatially resolved analysis which al- 
lows distinguishing between different construction times, 
functional areas, etc. A corresponding detailed study has 
been beyond the scope of this comparative investigation 
on different road networks, but is needed in order to bet- 
ter distinguish between the two modes ("natural" and 
"planned") of historical city growth. In particular, the 
considered sample of cities needs to be systematically ex- 
tended in order to gain substantial information about the 
influences of regional, topographic, and functional aspects 
on the properties of urban road networks. Corresponding 
studies will probably benefit from an additional consider- 
ation of the statistical properties of the associated dual 
graphs, which have been suggested to allow better distin- 
guishing between historical and planned urban areas. 

Finally, related to our observation of rather different 
large-scale morphologies of urban road systems, it will be 
relevant to study also network properties on larger scales, 
such as their modular structure, and relate the corre- 
sponding results to our findings presented in this work. 
One promising approach is utilising other measures of 
centrality than the node degree, such as node and link 
betweenness, which provide alternative characterisations 
of the importance of individual nodes and links in the 
network. We emphasise that unlike for the simple geomet- 
ric properties studied here, it is likely that the degree of 
similarity (or even universality) in the frequency distribu- 
tions of such higher-order network properties is consider- 
ably lower. A detailed corresponding investigation will be 
subject of future work. 
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